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We measure the time-dependent CP asymmetry parameters in B° — > K + K~ K° based on a 
data sample of approximately 227 million B-meson pairs recorded at the Y(4S) resonance with the 
BABAR detector at the PEP-II _B-meson Factory at SLAC. We reconstruct two-body B° decays to 
(j>{W20)K^ and 0(1020) K\, and the three-body decay K+ K~ K° s with <j>{lQ2Q)K° s excluded. For 
the B° -» (j>K° decays, we measure sin 2 (3 cB ((j)K ) = +0.50 ± 0.25 (stat)±°;°£(syst). The B° -» 
K + K~ K% decays are dominated by K + K~ S-wave, as determined from an angular analysis; we 
measure sin 2f3 cB (K + K~ K° s ) = +0.55 ± 0.22 (stat) ± 0.04 (syst) ±0.11 (CP), where the last error is 
due to the uncertainty in the fraction of CP-even contributions to the decay amplitude. We find no 
evidence for direct CP violation. 

PACS numbers: 13.25.Hw, 12.15.Hh, 11.30.Er 



In the Standard Model (SM) of particle physics, the de- 
cays B° -> K+K-K Q are dominated by b — > sss glu- 
onic penguin amplitudes, but can also be affected by am- 
plitudes that are suppressed by elements of the Cabibbo- 
Kobayashi-Maskawa (CKM) quark mixing matrix 
These CKM-suppressed amplitudes cannot be precisely 
known in a model- independent way Q , but are in general 
expected to be small |4J. Let 2(3 e s be the CP- violating 
phase difference between decays with and without mix- 
ing, and f3 = arg (— V c dV* b /V t dV t * b ) where Vij are elements 
of the CKM quark mixing matrix. The difference \j3— /3 e ff| 
is expected to be nearly zero, with theoretical uncertain- 
ties of a few degrees for B° — > 4>K Q . Larger uncertain- 
ties exist for B° -> K+ K~ K% with B° -> (j>K° s decays 
excluded, due in part to an extra CKM-suppressed tree 
amplitude contribution p|. 



Since additional decay diagrams with non-SM particles 
and interactions introducing new CP-violating phases 
may contribute to e g, measurements of sin2/3 e ff in 
these channels and their comparisons with the SM ex- 
pectation are sensitive probes for physics beyond the 
SM The value of sin 2(3 has been measured in 

B° -> J/i/jK° H0 with an average of 0.742 + 0.037. The 
BABAR and Belle collaborations have measured sin 2/3-2 
in (j)K° (+0.47 ± 0.34+°;°^ with 114 million BB pairs 



(B°) is given by 



and -0.96 ± 0.50 



-0.06 _ 

_ n with 152 million BB pairs((f>Kg 



+0.09 



only) respectively), and in K + K K® excluding 4>K^ 
(+0.57 ± 0.26 ± 0.04±g;^ with 122 million BB pairs 
and 

respectively) 



o.oo 

0.51±Q.26±0.05lg;^ with 152 million BB pairs 



At B factories, the neutral B mesons are exclusively 
produced in pairs. We select events for which one B 
(B rec ) is reconstructed as B° — > K + K~K° and the other 
(•Stag) is partially reconstructed as either B° or B°. We 
define At = t ICC — ttag to be the difference between the 
proper decay times of the B mesons. The decay rate 
f+ (f_ ) for the final state / when the Stag decays as a B° 



f±(At) 



-|At|/T B0 
4r R o 



[1± 



S f shx(Am d At) + C f cos(Am d At)}, (1) 



where r^o is the B° lifetime and Arrid is the B°-B° mix- 
ing frequency. The parameter Sf is non-zero if there 
is CP violation in the interference between decays with 
and without mixing, while a non-zero value for C / would 
entail direct CP violation. In the limit where the CKM- 
suppressed amplitudes do not contribute, the SM pre- 
dicts no direct CP violation (Cf = 0) since the domi- 
nant decay amplitudes have the same CP-violating phase, 
and that Sf = -r) f x sin2/3 off . For B° -> (f>K° de- 
cays, the effective eigenvalue i]f = —1, for B° — > 4>K® 
T) f = +1. For B° -» K+K-K° s decays, Vf = 2/ ovcn - 1, 
where / cvcn is the fraction of CP-even contributions to 
the B° — * K + K~K® amplitude. Then the value of T)t 
depends on the angular momentum of the K + K~ sys- 
tem: it is -1 for relative P-wave and +1 for S'-wave. 

In this paper, we present a measurement of sin 2(3 c s 
with almost twice the number of events as for the pre- 
vious BABAR results |8|, LlJ ■ We reconstruct B° candi- 
dates in two independent modes, <j)K° (with the K° ei- 
ther a K° or a Kfj and K + K~ K Q S (with the cp mass 
region excluded). K®'s are detected via their n + ir~ de- 
cay only. We extract the CP asymmetry parameters us- 
ing extended maximum-likelihood fits. Using an angu- 
lar moment analysis [ll|, we extract the K + K~ P-wave 
fractions in the data. These fractions are used to check 
the assumption that 77/ = —1 for <pK2, and +1 for <j)K® 
by bounding the 5-wave contamination in the <f> mass 
region, and to measure r\f for K + K~ K®. 

This analysis is based on 227 million BB pairs collected 
with the BABAR detector [l2J at the PEP-II asymmetric- 
energy e + e~ storage rings at SLAC, operating at the 
T(4S) resonance (center-of-mass (cm.) energy y/s = 
10.58 GeV). In Ref. we describe the silicon vertex 
tracker (SVT) and drift chamber (DCH) used for track 
and vertex reconstruction, the electromagnetic calorimc- 
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ter (EMC) and instrumented flux return (IFR) used for 
K® reconstruction, and the detector of internally re- 
flected Cherenkov light (DIRC), which, together with the 
EMC, the IFR, and the ionization dE/dx from the SVT 
and DCH, is used for particle identification. 

The -B°-candidate reconstruction and selection is sim- 
ilar to that described in Refs. and We consider a 
K + K~ pair to be a <f> candidate if its invariant mass is 
within 15MeV/c 2 (about three times the apparent width 
in the K + K~ invariant mass spectrum) of the central 
4> mass value 13] . For a given B°B° meson pair, we 
obtain At from the measured distance between the fully 
reconstructed B rcc meson decay point and the B ta , g decay 
point along the beam direction, and the known boost of 
the T(45) system (/?7=0.56). A multivariate tagging al- 
gorithm determines the flavor of the -Btag meson [6j and 
classifies it in one of seven mutually exclusive tagging 
categories. 

We use two kinematic variables to discriminate be- 
tween signal B decays and combinatorial background. 
The energy difference between the measured e + e~ cm. 
energy of the B candidate and y/s/2 is AE. Its dis- 
tribution peaks at zero for signal, with a width of 
about 20MeV for (f>K° and K+K~K°. The width 
is only about 3MeV for <pK®, because for this mode 
we constrain the B Q candidate's mass to the nominal 
value 13]. The beam-energy-substituted mass, toes, 
is used for candidates without a K®. It is defined as 
toes = \/(s/2 + Pi -PB) 2 /E'f -p%, where the B mo- 
mentum p b and the four-momentum of the initial state 
(Ei, pi) arc defined in the laboratory frame. It peaks 
at the B° mass for signal, with a width of about 3 MeV. 
For (f>K® candidates, we require \AE\ < 100 MeV and 
toes > 5.21 GeV/c 2 ; for 4>K® candidates, we require 
| AS | < 80 MeV; and for K + K~K° S candidates, we re- 
quire \AE\ < 200 MeV and m ES > 5.2 GeV/c 2 . 

The dominant background is continuum e + e~ — > qq 
(q = u, d, s, c) events; these tend to be jet-like in the 
e + e~ cm. frame, while B decays tend to be spherical. 
To enhance discrimination between signal and contin- 
uum, we use Fisher discriminants (T) to combine four 
event-shape-related variables 0, The other back- 

ground originates from B decays. For the 4>K° final 
state, opposite-CP contributions from the K + K~K° fi- 
nal state (K + K~ S-w&ve) are estimated from data with 
a moment analysis ^lj (see below) to be less than 6.6% 
at 95% confidence level. The mode <j>K® has additional 
background. Its dominant CP contamination is from the 
mode (f>K*° — > tfiK®ir°, for which we expect approxi- 
mately eight events in the region \AE\ < 10 MeV. In 
the final likelihood fit we explicitly parameterize back- 
grounds from B decays both with and without charm. 

For the K + K~ K® mode, we apply invariant mass 
cuts to suppress background from B decays that proceed 
through a b — > c transition, namely those containing D , 
J/tp, Xco, or tp(2S) decaying into K + K~, or D + or Df 



decaying into K + K®. Finally, to suppress B decays into 
final states with pions, we require the rate for a charged 
pion to be misidentified as a kaon to be less than 2%. 

A total of 4,300 , 8,238, and 27,368 events have a can- 
didate that passes the 4>K®, 4>K®, or K + K~K® selection 
criteria, respectively. From simulation, we find the final 
selection efficiencies for signal to be 40%, 20% and 26%, 
respectively. 

We extract the K + K~K°, <j>K°, and cj>K° event 
yields and CP parameters with two extended maximum- 
likelihood fits. One is to the K + K~K® candidates; the 
other is to both the <j>K® and <j>K® candidates, with the 
assumption that C^o — C<t>K\ and = — . We 

verified the fit procedure for the 4>K° mode with samples 
of 4>K + and J/ipK® events. We found for the former a 
null asymmetry as expected, and for the latter results 
that are consistent with previous measurements 0. We 
verified the fit procedure for the K + K~K® mode with 
a sample of K®K®K + events, for which we found a null 
asymmetry as expected. 

The likelihood function used in each extended 
maximum-likelihood fit to its Nk candidates tagged in 
category k is 

N k , lis \ 

i=l j=l ' 

(2) 

where N' k is the sum of the signal, continuum, and ns B- 
background yields tagged in category k; Ns is the num- 
ber of 4>K®, <fiK®, or K + K~K® signal events; efc is the 
fraction of signal events tagged in category k; Nc,k is the 
number of continuum background events tagged in cat- 
egory k; Nbj is the number of S-background events of 
class j; and e^fc is the fraction of B-background events 
of class j tagged in category k. Each i3-background class 
comprises similar B decays. The B-background event 
yields are fixed parameters and are zero for the <f>Kg sam- 
ple. The total likelihood L is the product of the likeli- 
hoods for each tagging category. 

The probability density functions (PDFs) , V% , 
and V^ k , for signal, continuum background, and B- 
background class j, respectively, are the products of the 
PDFs of the discriminating variables. The signal PDF 
is thus given for the K + K~K® sample by 'P(toes) • 
V(AE) ■ V{T) ■ V(At;a At ), for the (f>K° s sample by 
V{m ES ) -V(AE) -V{T) -V{m KK ) -V{cos 9 H )-V{At; a At ), 
and for the <j)K° L sample by V{AE) ■ V{T) ■ V(m KK ) ■ 
■p(cos 6h) 'V{At; (TAt), where Or is the angle between the 
K + candidate and the parent i? roc flight direction in the 
K + K~ rest frame. The quantity a At is the uncertainty 
in the measurement of At for a given event. The time- 
dependent CP parameters defined in Eq. JU, diluted by 
the effects of mistagging and the At resolution, are con- 
tained in {At, a At)- As in our J/tpK® analysis [||, 
the At-resolution function for signal and _B-background 
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events is a sum of three Gaussian distributions, which 
have two distinct means as well as three distinct widths. 
The widths are the error of the measured At scaled by 
three independent factors. 

In the fits to data, we leave unconstrained the param- 
eters describing the CP asymmetry, the Ai-resolution 
functions, the tagging characteristics, and the event 
yields. We also leave unconstrained the means of the sig- 
nal 77iES and AE Gaussian PDFs, the widths of the sig- 
nal AE PDFs, the mean of the signal tukk PDF (which 
is parameterized by a relativistic P-wave Breit-Wigner 
function), and all parameters of the K + K~K® candi- 
dates' signal PDF for T. We take from simulation any 
other parameters of the toes, AE, T, cos 9h, and tukk 
PDFs for signal and B background. The parameters de- 
scribing the signal and P-background At-resolution func- 
tion are determined by a simultaneous fit to an indepen- 
dent sample of reconstructed B° decays to flavor eigen- 
states, with more than 100,000 events M. We use the 
world-averaged values for t b o and Ama [13j . The fits to 
the (j>K° and K + K~K2, candidates have a total of 35 and 
34 free parameters, respectively. 

We use an angular moment analysis based on the 
cosOh distribution to extract the K + K~K® CP content, 
and also to bound the 5-wave contamination in the <f) 
mass region. In this approach, we expand the decay dis- 
tribution for a given K + K~ invariant mass in terms of 
moments (Pi) of conveniently normalized Legendre poly- 
nomials Pe(cos8n)- 

\A(m KK )\ 2 = J2( P i)- p dcose H ), (3) 

i 

where A(m,KK) is the mass-dependent decay ampli- 
tude. We normalize P^(cos6*#) such that the integral 
of Pf(cos 9h) 2 over cos9h from —1 to 1 equals unity. We 
extract the moments by summing over all events: 

(p t ) = ^/n<-u^//jiv, ,. (4) 

3 

where Wj is the weight for event j to belong to the si gnal 
decay and is calculated by the sPlot technique of Ref . [lj. 
The efficiency Ej is evaluated from a large MC sample 
in bins of ttlkk and cos 6h ■ Limiting ourselves to the 
two lowest partial waves, we can write the total decay 
amplitude in terms of the S-wave (CP-even) and the P- 
wave (CP-odd) amplitudes, 

A(m K K) ~ As{m KK )P o (cos0 H ) + 

e**'Ap(m K K)Pi(coB6 H ), (5) 

where (f> p is the relative phase between the real partial- 
wave amplitudes As(ttikk) and Ap(vtikk)- If we com- 
pare Eq. JSJ to Eq. ©, we can relate the moments (of 
order i < 2) to the wave intensities and thus to the total 



fraction of CP-even events, / cven , as 

, _ A s (m KK ) 2 

As{mKK) 2 + Ap(rriKK) 2 

= 1 ~vw (6) 

where As(ni>KK ) 2 and Ap(rriKK ) 2 are the S- and P-wave 
intensities, respectively. In the normalization, the total 
number of signal events is \/2(Pq). 

TABLE I: Systematic uncertainties on the CP parameters. 



Source 






Skkk 


Ckkk 


Detector effects 


±0.02 


±0.02 


±0.02 


±0.01 


DCSD 


±0.01 


±0.03 


±0.00 


±0.03 


Fit bias 


±0.01 


±0.01 


±0.02 


±0.01 


B°-B° tagging 


±0.01 


±0.02 


±0.00 


±0.01 


.S-wave contamination +0.06 ±0.02 






Other 


±0.03 


±0.02 


±0.01 


±0.01 


Total 


+U.U7 
-0.04 


±0.05 


±0.03 


±0.04 



Systematic errors on the CP-asymmetry parameters 
are listed in Table |I] We account for uncertainties in 
the At resolution, the beam-spot position, and the de- 
tector alignment. We also estimate errors due to the 
effect of doubly CKM-suppressed decays (DCSD) of the 
^tag 0|- The uncertainty due to possible biases in the 
fit procedure is conservative and includes effects on the 
CP parameters of correlations among the fit variables, 
which have been determined with full-detector MC sim- 
ulations. Uncertainties in the B°-B° tagging efficiency in 
both signal and background are also included. Finally, 
we account for errors due to the CP content of the back- 
ground, uncertainties in the PDF parameterization, and 
the uncertainties of t b o and Ama 13j . For each mode we 
add the individual contributions in quadrature to obtain 
the total systematic uncertainty. 

We also consider the systematic error due to the 
CP-even fraction of the K + K~K® mode. We do not 
find evidence for the existence of higher moments (Pi), 
£ = 3 ... 6, which could arise from intermediate P-wave 
decays into K + K~ or decays proceeding through an 
isospin-1 resonance into K K%. Nevertheless, we esti- 
mate a systematic error from the P-wave by examining 
the (P2) moment in the K + K~ mass region (1.1—1.7) 
GeV/c 2 , corresponding to the / 2 (1270), a 2 (1320)°, and 
/2(1525) resonances, and assuming that (P2) arises only 
from P-wave and S-D interference. Since the moment it- 
self is consistent with zero, we assign a systematic error of 
4% based on the (P2) error. We account for the possible 
presence of a (980)+, a (1450)+, and a 2 (1320)+ in the 
K ± Kg subsystem (4.6%). We also estimate a bias due 
to the modeling of the efficiency from MC events (2.5%). 
We find the total systematic error on / OV cn to be ±0.06. 
This leads to a systematic error on sin2 / 9 e gf of ±0.11. 
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Table[n]shows the measured CP parameters and yields 
from the final extended maximum- likelihood fits. Note 
that when fitting sin2/3 e fj for K + K~ K®, we constrain 
Ckkk to zero. All yields are consistent with our pre- 
viously measured branching fractions 0, 0] . Figure Q] 
shows the signal-enhanced distributions of toes for <f>K® 
and K + K~K® events and of AE for (f>K® events, to- 
gether with the result from the final extended maximum- 
likelihood fits. Figure shows the time-dependent asym- 
metry distributions. As a cross check, we also fit <fiK® 
and 4>K° L separately. Our fit to only <j>K® events gives 
S = 0.29 ± 0.31 and C = -0.07 ± 0.27. Our fit to only 
cj)K° events gives S = 1.05 ± 0.51 and C = 0.31 ± 0.49. 



c 

CD 
> 
LU 



O 

> 

CD 
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5.26 
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FIG. 1: Distributions of (a) jties for 4>Ks candidates, (b) AE 
for <\>K% candidates, and (c) iti-es for K + K~ K® candidates 
excluding 4>K%, together with the results from the final ex- 
tended maximum-likelihood fits after applying a requirement 
on the ratio of signal likelihood to signal- plus-background like- 
lihood (computed without the displayed variable) to reduce 
the background. The requirement is chosen to roughly max- 
imize Ns/(Ns + Nc) where Nc is the total number of con- 
tinuum events, and is applied only for the purpose of making 
these plots. The curves are projections from the likelihood 
fits for total yield (solid lines), continuum background (short 
dashed lines), and total background (long dashes in (b) only). 
The efficiency of the likelihood-ratio cut is (a) 79% for sig- 
nal and 5% for background, (b) 35% for signal, 16% for B- 
background, and 3% for continuum background, and (c) 77% 
for signal and 5% for background. 

For the K + K~K® final state including the <f> mass re- 
gion, the distributions of the S- and P-wave intensities, 
and the CP-even fraction, as a function of K + K~ in- 
variant mass, are shown in Fig. |3J The total fraction 
of CP-even events with the <f> mass region excluded is 



TABhE II: CP-asymmetry parameters and yields from the 
final extended maximum-likelihood fits, as well as the fraction 
of CP-even contributions to the amplitude, / evcll , which is 
assumed to be zero for 4>K® and and unity cj>K°. The first 
errors are statistical, and the second are systematic; the third 
error on sin2/3 c ff for K + K~ K® is due to the uncertainty in 
the CP content. The values of S and C are fit simultaneously 
for the 4>Ks an d <t>K° candidates; the sign of S for 4>K% is 
shown. When finding sin2/3 c ff for K + K~K%, we constrain 
Ckkk to 0. 





cAA'" 


K + K~K U S 




<f>K° s 4>K° L 


(no <t>K°) 


sin 2(3 C B 


+0.50 ± 0.25t££ 


+0.55 + 0.22 + 0.04 + 0.11 


/even 


1 


0.89 + 0.08 + 0.06 


S 


+0.50 ± 0.25±g"gl 


-0.42 + 0.17 + 0.03 


C 


0.00 + 0.23 + 0.05 


+0.10 + 0.14 + 0.04 


Yield 


114 + 12 98 + 18 


452 ± 28 



1 

0.5 


-0.5 
-1 
1 

0.5 


-0.5 
-1 

0.5 


-0.5 






-6 -4 -2 2 4 6 
At ( P S) 

FIG. 2: The time-dependent asymmetry distributions for (a) 
4>K%, (b) cj>K L , and (c) K + K~ K° s with no <f>K° s decays. The 
asymmetry is defined as A g0 /go = {N B o —Ng ) / {N B o +Ng ) , 
where N B o (Ngo ) is the number of Stag mesons identified as 
a B°(B°) for a given measured value of At. The signal-to- 
background ratio is enhanced with a cut on the likelihood 
ratio as in Fig. Q 

given in Table [H] We successfully verified our value of 
/even with a different method [17j that uses the event 
rates in B + — > K + K®K® and the isospin-related channel 
B° -> K+K-K°. 

To summarize, in a sample of 227 million BB meson 
pairs, we measure the CP content and CP parameters 
in P°-meson decays into <pK°, and into K + K~K® with 
the <j) mass region excluded. We determine the fraction 
of CP-even and CP-odd contributions with an angular 
analysis. In B° — > <pK°, our values for sin2/3 e ff and C^k 
are in good agreement with our previously published val- 
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FIG. 3: Distributions of S- and P-wave intensities and CP 
even fraction as a function of K + K~ invariant mass. Notice 
that the first bin integrates a wider mass range than the <j> 
resonance occupies. Insets show S- and P-wave intensities in 
the <j> mass region. 

ues and the small S'-wave contamination is treated 
as a systematic uncertainty. In B° — > K + K~K®, the 
K + K~ system is observed to be dominated by S'-wave; 
this, along with the measured value of sin2/3 e ff> is consis- 
tent with previous measurements based on isospin sym- 
metry |l0| . Both of our sin 2f3 c $ values are consistent 
to within one standard deviation with the value of sin 2(3 
measured in £>° — > ccs decays Q. 
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